Shape anisotropy of polymers in disordered environment 
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We study the influence of structural obstacles in a disordered environment on the size and shape 
characteristics of long flexible polymer macromolecules. We use the model of self-avoiding random 
walks on diluted regular lattices at the percolation threshold in space dimensions d = 2, 3. Applying 
the Pruned-Enriched Rosenbluth Method (PERM), we numerically estimate rotationally invariant 
universal quantities such as the averaged asphericity {Ad) and prolateness {S) of polymer chain 
configurations. Our results quantitatively reveal the extent of anisotropy of macromolecules due to 
<N ' the presence of structural defects. 
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I. INTRODUCTION 

*?! 

Topological properties of macromolecules, such as the shape and size of a typical polymer chain configuration, are 
^ of interest in various respects. The shape of proteins affects their folding dynamics and motion in a cell and is relevant 
in comprehending complex cellular phenomena, such as catalytic activity ll[ . The hydrodynamics of polymer fluids 
is essentially affected by the size and shape of individual macromolecules [2j; the polymer shape plays an important 
role in determining its molecular weight in gel filtration chromatography [3] . 

Already in 1934 it was realized [4| that viscous properties of polymer solutions are significantly different than 
predicted by the theory for dissolved sphere-like molecules: flexible polymer chains in good solvents form crumpled 
coil shapes, which are surprisingly anisotropic. Since then, a considerable amount of work has been done in exploring 
size and shape characteristics of various macromolecules |5l-ll8l| . 

Sole and Stockmayer Q introduced as a shape measure of macromolecules the normalized average eigenvalues 
£N| . Xi of the gyration tensor. Numerical simulations in d = 3 dimensions give {(Ai), (A 2 ), (A 3 )}={0. 790, 0.161, 0.054} 
[H, indicating a high anisotropy of typical polymer configurations compared with the purely isotropic case {1/3, 
1/3, 1/3}. To compute the quantities Ai analytically is, however, difficult, because one must explicitly diagonalize 
the gyration tensor for each realization in an ensemble of polymers. It was proposed [ToL [Tlj to characterize the 
asymmetry of polymer configurations by rotationally invariant universal quantities, such as the averaged asphericity 
• ■ (Ad) and prolateness (S). {Ad) takes on a maximum value of one for a completely stretched, rod-like configuration, 
and equals zero for spherical form, thus obeying the inequality: < {Ad) < 1. The quantity (S), defined in d = 3, 
takes on a positive value for prolate ellipsoid-like configurations, and is negative for oblate shapes, being bounded to 
the interval — 1/4 < (S) < 2. To characterize the size measure of a single flexible polymer chain, one usually considers 
the mean-squared end-to-end distance (R 2 ) and radius of gyration (Rq), both governed by the same scaling law: 
(i? 2 ) ~ {R%) ~ N 2l/ , where N is the mass of the macromolecule (number of monomers in a polymer chain) and v a 
universal exponent (y > 1/2 (d < 4), v — 1/2 [d > 4)). The ratio of these two characteristic distances, the so-called 
size ratio Qd = (R 2 e )/(R G ), also appears to be a universal, rotationally invariant quantity {gd > 6 (d < 4), gd = 6 

(d>4))0- 

Numerous studies indicate that a typical flexible polymer chain in good solvent takes on the shape of an elongated, 
prolate ellipsoid. In particular, using the data of x-ray crystallography and cryo-electron microscopy, it was found 
that the majority of nonglobular proteins are characterized by A3 values from 0.5 to 0.7 and S values from to 0.6 
T7I ill. The shape parameters of polymers were analyzed analytically, based on the renormalization group approach 
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9l. ll0lll6| . and estimated in numerical simulations [El H EH [HI- Previous estimates of the shape and size characteristics 
of flexible polymer chains in d = 2, 3 are given in Table [2 

In real physical processes, one is often interested in the question as to how structural obstacles (impurities) in the 
environment alter the behavior of a system. The density fluctuations of obstacles lead to a large spatial inhomo- 
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geneity and create pore spaces, which are often of fractal structure [l9j . In polymer physics, of great importance 
is understanding of the behavior of macromolecules in the presence of structural disorder, e.g., in colloidal solutions 
20] or microporous membranes [2lJ. In p articular, a related problem is relevant when studying the protein folding 
dynamics in the cellular environment [22|. Biological cells can be described as a highly disordered environment due 
to the presence of a large amount of soluble and insoluble biochemical species, which occupy up to 40% of the total 
aquabased volume Structural obstacles strongly effect the protein folding and aggregation [24T - I271 ] . Recently, 

it was realized experimentally [28j that macromolecular crowding has a dramatic effect on the shape properties of 
proteins. To explain the physics behind the macromolecular crowding effects, a statistical theory of excluded volume 
interaction between proteins and hard particles was developed j29( . The folding of proteins in crowded environment 
was studied numerically by off-lattice polypeptide chain simulations in the presence of repulsive spherical particles in 

Refs. QmJ cm 

In the language of lattice models, the disordered environment with structural obstacles can be considered as a 
lattice where some amount of randomly chosen sites contain defects, which are to be avoided by the polymer chain. 
Of particular interest is the case, when the concentration of lattice sites allowed for the polymer chain equals the 
critical concentration p c — 0.592746 (d — 2) [32], p c = 0.31160 (d — 3) [33[ and the lattice becomes percolative: a 
percolation cluster, having a fractal structure, occurs 34]. The lattice model with percolation has played a key role 
in statistical physics for describing structurally disordered systems. This approach is limited, however, in the realm 
of biological systems and processes. The core of this difficulty is the highly specialized and nonequilibrium nature 
of the biological world. Still, some aspects such as the universal configurational properties of bioproteins seem to be 
amenable to analyses via this simplified model [22|, [26[ . 

Studying processes on percolative lattices, one encounters two possibilities. In the first, one considers only per- 
colation clusters with linear size much larger than the typical length of the physical phenomenon under discussion 
(polymer chain length in our case). The other statistical ensemble includes all lattice sites free of defects, which can 
be found in a percolative lattice. In the latter case, the polymer can be trapped (localized) in confined regions of 
pure sites, and so-called localization phenomena occur which lead to decreasing the size of the macromolecule. Such a 
situation has been studied in Ref. [131 ] , realizing the shrinking of polymer shapes with a trend to decreasing anisotropy. 
In what follows, we will be interested in the former case, when a polymer chain resides only on the percolation cluster 
of fractal structure; this results in increasing the swelling of the polymer coil compared with the pure solution case 
[35M40l ] . However, the important question of how do the shape parameters of a polymer chain change quantitatively, 
when the polymer is located on a fractal cluster, still is completely unresolved. 

The purpose of the present paper is to report computer simulations of a model of flexible polymer chains on the 
backbone of fractal percolation clusters in d = 2, 3. We aim to obtain numerical estimates of (A^ c ), (S Po ) and g^", 
and thus quantitatively describe the change in asymmetry of typical polymer configurations due to the presence of 
structural obstacles. The rest of the paper is organized as follows: in the next section we introduce the discretized 
mathematical model of polymer chain and Section III describes the details of our computer simulations. We give 
discussions of our results in Section IV and end up by giving conclusions and an outlook. 



II. DESCRIPTION OF THE MODEL 



Let R n = {x T \, . . . , x^} be the position vector of the nth monomer of a polymer chain (n = 1. . . . , N). The measure 
of the shape properties of a specified spatial conformation of the chain can be characterized @, GU El m terms of the 
gyration tensor Q with components: 

1 N 

Qij — Jj£_j^ X%n ~ X CM)( X n ~ x Cm)i *'i = l'---'^' (1) 

n=l 

with x l CM — XmLi X n/N being the coordinates of the center-of-mass position vector Rcm- 

The spread in eigenvalues of the gyration tensor describes the distribution of monomers inside the polymer 
coil and thus measures the asymmetry of a molecule; in particular, for a symmetric (spherical) configuration all the 



TABLE I: Size ratio, averaged asphericity and prolateness of flexible polymer chains on regular lattices, a: Ref. [jj, 6: Ref. 
[12, c: Ref. [3. 

d ffd (A d ) (S) 

2 7.14±0.03 a 0.501 ± O.OOS 6 

3 6.249 ± 0.03 c 0.431 ± 0.002 c 0.541 ± 0.004 c 
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eigenvalues A^ are equal. 

It was found convenient to characterize the shape of polymers by rotationally invariant universal combinations of 
components of the gyration tensor [HI EH ■ The first invariant of Q is the squared radius of gyration 

^ N d 

R 2 G = -J2^~ ^M? = J2 Q« = Tr Q> ( 2 ) 

n—1 i—1 

which measures the average distribution of monomers with respect to the center of mass. Let A = TrQ/d be 
the average eigenvalue of the gyration tensor. Then the extent of asphericity of a polymer chain configuration is 
characterized by the quantity Ad defined as fiol ]: 

1 >f (A t -A) 2 d TrQ 2 
d d{d-l){r[ A 2 d-l(TrQ)2< W 

with Q = Q — A I (here I is the unity matrix). This universal quantity equals zero for a spherical configuration, where 
all the eigenvalues are equal, Xi — A, and takes a maximum value of one in the case of a rod-like configuration, where 
all the eigenvalues equal zero except of one. Thus, the inequality holds: < Ad < 1. Another rotationally invariant 
quantity, defined in three dimensions, is the so-called prolateness S fiol. ITij: 

q Ilti(A.-A) _ detQ 

If the polymer is absolutely prolate, rod-like (Ai ^ 0, A2 = A3 = 0), it is easy to see that S equals two. For absolutely 
oblate, disk-like conformations (Ai = A2, A3 — 0), this quantity takes on a value of —1/4. In general, S is positive for 
prolate ellipsoid-like polymer conformations (Ai >• A 2 ~ A3) and negative for oblate ones (Ai «A 2 » A3), whereas its 
magnitude measures how oblate or prolate the polymer is. Note that since A and the quantities in ©-O are expressed 
in terms of rotationally invariants, there is no need to explicitly determine the eigenvalues A^ which simplifies the 
numerics significantly. 

The average of quantities ©-Q for a given polymer chain length N, denoted as (...), is performed over an ensemble 
of possible configurations of a chain. Note that some analytical and numerical treatments avoid the averaging of the 
ratio in ([3]), (0| and evaluate quantities: 

which should be distinguished from the averaged asphericity and prolateness: 




nti(A .-A) 

d(d-l)\U A 2 /' \ A 3 



(A d ) = -^-AY:-^f^), (S) = ( ^_: ^ ■ (6) 



Contrary to {Ad) and (S), the quantities ([5]) have no direct relation to the probability distribution of the shape 



parameters Ad and S. As pointed out by Cannon et al. 15|, this definition overestimates the influence of larger 
polymer configurations on the mean shape properties and suppresses the influence of compact ones. This artificially 
leads to overestimated values for shape parameters. The difference between (Ad) and Ad on regular lattices was found 
to be really large (A 2 = 0.625 ± 0.008, A 3 = 0.546 ± 0.008 [2], which should be compared with the data in Table H}. 



III. THE METHOD 



We start with regular lattices with sites assigned to contain an obstacle with probability 1— p c and be allowed for the 
polymer chain otherwise. To obtain the backbone of a percolation cluster on a given disordered lattice, we apply an 
algorithm consisting of the following two steps: first finding the percolation cluster based on the site-labeling method 
of Hoshen and Kopelman [iilj and then extracting the backbone of this cluster [42j (the algorithm is explained in detail 
in our previous papers (40j). Since the self-avoiding walk trajectory can be trapped in dangling ends of percolation 
clusters, the infinitely long chains can only exist on the backbone of the cluster. The question about differences in 
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universal configurational properties of SAWs walking on a percolation cluster and its backbone was addressed recently, 
e.g., in Ref. [40] ■ We constructed 1000 clusters in each space dimension. 

To study shape properties of typical polymer chain configurations modelled by self-avoiding random walks (SAWs) 
on the constructed percolation clusters, we use the pruned-enriched Rosenbluth method (PERM 1 ) |43lj. combining 
the original Rosenbluth- Rosenbluth (RR) algorithm of growing chains [4J| and population control [45[ . The growth 
process starts at the center of the percolation cluster, and each nth monomer is placed at a randomly chosen neighbor 
site of the last placed (n — l)th monomer (n < N, where N is total length of the polymer). If this randomly chosen site 
is already visited by a chain trajectory or does not belong to the percolation cluster, it is avoided without discarding 
the chain, but the bias is corrected by means of giving a weight W n ~ (IIzLi TO t° each sample configuration at 
the nth step, where m; is the number of free lattice sites to place the Zth monomer. The growth is stopped when 
the total length N of the chain is reached, then the next chain is started to grow from the same starting point. The 
configurational averaging for any quantity of interest then has the form: 



((•■•)) 



conf 
N • 



(7) 



The Rosenbluth method, however, also suffers from attrition: if all next neighbors at some step (n < N) are 
occupied, i.e., the chain is running into a "dead end", the complete chain has to be discarded and the growth 
process has to be restarted. Grassberger [43j proposed a considerable improvement of the efficiency by increasing the 
number of successfully generated chains. The weight fluctuations of the growing chain are suppressed in PERM by 
pruning configurations with too small weights, and by enriching the sample with copies of high- weight configurations. 
These copies are made while the chain is growing, and continue to grow independently of each other. Pruning and 
enrichment are performed by choosing thresholds and depending on the estimate of the partition sum for 
the n-monomer chain. These thresholds are continuously updated as the simulation progresses. If the current weight 
W n of an n-monomer chain is less than W<, the chain is discarded with probability 1/2, otherwise it is kept and its 
weight is doubled. If W n exceeds W> , the configuration is doubled and the weight of each identical copy is taken 
as half the original weight. For a value of the weight lying between the thresholds, the chain is simply continued 
without enriching or pruning the sample. For updating the threshold values we apply similar rules as in |46l |47| : 
= C(Z n /Z 1 )(c n /ci) 2 and W 7 f — 0.2W>, where c„ denotes the number of created chains having length n, and 
the parameter C controls the pruning-enrichment statistics. After a certain number of chains of total length N is 
produced, the given tour is finished and a new one starts. The pru ning-enrichment control is adjusted such that on 
average 10 chains of total length N are generated per each tour [47j. 

For estimations of quantities of interest we have to perform two types of averaging: the first over all polymer con- 
figurations on a single percolation cluster according to ([7]); the second average is carried out over different realizations 
of disorder, i.e., over all percolation clusters constructed: 



1 M 



(8) 



where M is the number of different clusters and the subscript i means that a given quantity is calculated on the 
cluster i. 
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FIG. 1: Correlation between A3 and S values (a) and between Rq and R\ values (b) in d = 3 for chain length TV = 120. 
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FIG. 2: Averaged asphericity of polymer configurations in d — 2 (a), d — 3 (b) and averaged prolateness in d = 3 (c). Lower 
lines: pure lattice, upper lines: percolation cluster. Dashed lines show the results of least-square fitting with the ansatz 



IV. RESULTS 



We construct percolative lattices of edge lengths up to L max =400, 200 in dimensions d=2, 3, respectively, and 
estimate the mean shape parameters ©• The disorder averaging is performed over 1000 percolation clusters in each 
space dimension. 

At first, let us analyze the connections (correlations) between values of the shape characteristics Ad and 5*. In Fig. 
[IJa) we present our data for simultaneous estimates of these quantities for each N = 120-step SAW configuration 
of the ensemble in d = 3. As it is clear from the definitions ([3]) and (U]), for spherical configurations both A3 and 
S3 are close to zero. S equals zero also for the cases when Ai + A2 = 2 A3, corresponding to a non-zero asphericity 
value. Increasing positive values of S describe increasing the elongation and tending to rod-like structures, which 
corresponds also to increasing A3 values. The negative values of S describe the non-spherical oblate structures, also 
corresponding to non-zero A3 values. In Fig. []Tb) we give also our data for simultaneous estimates of R 2 G and R\ in 
an ensemble of N = 120-step SAWs in d — 3; these quantities are less nicely correlated. 

Figures [2] and [3] present simulation data for (Ad), (S) and gd = (RI)/(Rq) as functions of the chain length N in 
d = 2, 3. For comparison, we also evaluated the shape parameters on the pure lattices of the same edge lengths. Note 
that the size ratio gd is rather a delicate quantity due to large fluctuations in the sample. 

For the case of the pure lattice, it is evident that the shortest chains are very elongated in space and they become a 
little more spherical with increasing chain length; thus (Ad(N)) and (S(N)) are decreasing functions of chain length 
N. On a percolative lattice, these quantities behave in a different manner: they increase gradually with increasing N. 
The structure of the fractal percolation cluster makes the longer polymer chain configurations to be more and more 
prolate. The principal qualitative conclusions which we can derive from Figs. [2] and [3] is that the shape parameters 
of typical polymer configurations change significantly relative to the obstacle-free case; the shape tends to be more 
anisotropic, elongated due to the fractal structure of the lattice. 

For finite chain length N, the values of shape parameters differ from those for infinitely long chains. This finite-size 




FIG. 3: Size ratio gd = (R\) / (R%) of polymer configurations in d = 2 (a) and d = 3 (b). Lower lines: pure lattice, upper lines: 
percolation cluster. Dashed lines show the results of least-square fitting with the ansatz ([9]). 
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FIG. 4: Probability distributions of asphericity in d — 2 (a, b) and d = 3 (c, d) and prolateness in d = 3 (e, /) of SAWs on a 
pure lattice (left column) and on a percolation cluster (right column) for chain length N = 120. 



deviation obeys scaling behavior with TV: 

(A d (N))^(A d )+b 1 N- A , 

(S(N)) = (S)+b 2 N- A , (9) 
g d (N)=g d + b 3 N~ A , 

where bi, b 2 , &3 are constants and A is the correction-to-scaling exponent: A(d = 2) = 1.5 [48j], A(d = 3) = 1.7 [49| , 
The shape parameter estimates can be obtained by least-square fitting of ([9]). For the case of the pure lattice, we receive 
a nice agreement with the existing data of Table Q] (A 2 ) = 0.506 ± 0.002, (A3) = 0.435 ± 0.002, (S) = 0.545 ± 0.002. 
Our results for the shape characteristics of SAWs on a percolation cluster are given in Table HH Note, that since we 
can construct percolative lattice only up to a finite size L, it is not possible to perform very long SAWs on it. For each 
L, the correct statistics holds only up to some "marginal" number of SAWs steps iV mar g ~ L 1 / 1/ saw ) with i*s c AW being 
the size scaling exponent on the percolation cluster [401 ] . We take this into account when analyzing the data obtained; 
for each lattice size we are interested only in values of N < -/V marg , thus avoiding distortions, caused by finite-lattice 
effects. We have also estimated the shape parameters defined by ratios of averages |0 with (. . .) replaced by (. . .) in 
the disordered case; as expected, the obtained values are considerably larger than averages of ratios, cf. Table ITT1 

The probability distributions of shape parameters in d — 2, 3 at fixed chain length N — 120 are given in Fig. |U 
The distribution functions for A 2 and A^" are rather unsymmetrical with the most probable value larger than the 
mean value. In d — 3, the asphericity distribution function is broad but quite symmetric, with the most probable 
and mean value being nearly equal. The distribution function for S is rather unsymmetrical with a most probable 
value of zero, which is shifted for S Pc . The shapes of these distributions indicate that the majority of polymer chain 
configurations has prolate asymmetry. 
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TABLE II: Our estimates for size ratio, averaged asphericity and prolateness of SAWs on percolation clusters. 



2 7.96 ± 0.01 0.571 ± 0.005 0.68 ± 0.01 



3 7.44 ±0.02 0.531 ±0.005 0.743 ±0.005 0.61 ± 0.01 0.97 ± 0.01 



V. 



CONCLUSIONS 



We study the universal size and shape characteristics of flexible polymer macromolecules in an environment with 
structural obstacles. The measure of the shape properties of a specified configuration of the polymer chain is char- 
acterized in terms of the gyration tensor Q. The rotationally invariant quantities, constructed as combinations of 
components of Q, such as the averaged asphericity (Ad) and prolateness (S) of typical chain realizations are of interest. 
Another quantity of interest is the universal size ratio gd = (R^)/(R%). We address the question, how the polymer 
shape anisotropy is quantitatively influenced by the presence of disorder, which is important in understanding many 
real physical processes. 

We use the lattice model of self-avoiding random walks on a disordered lattice exactly at the percolation threshold, 
when a percolation cluster with fractal structure emerges. Studying processes on percolative lattices, one encounters 
two possibilities. In the first, one considers only percolation clusters with linear size much larger than the typical 
length of the physical phenomenon under discussion (polymer chain length in our case). The other statistical ensemble 
includes all lattice sites free of defects, which can be found in a percolative lattice. We considered the former 
case, with a polymer chain residing only on the backbone of a percolation cluster. Applying the Pruned-Enriched 
Rosenbluth Method, we performed computer simulations in d = 2 and d = 3 and obtained numerical estimates for 
the averaged asphericity, prolateness and size ratio of self- avoiding walks on a percolation cluster. All the shape 
characteristics increase gradually with increasing polymer chain length; the structure of fractal percolation cluster 
makes the longer polymer chain configurations to be more and more prolate. Our results quantitatively indicate 
that the shape parameters of typical polymer configurations change significantly relative to the obstacle-free case; the 
shape tends to be more anisotropic and elongated due to the fractal structure of the lattice. 



Note that recently a related model has been studied analytically in Ref. [50( , where the shape properties of polymers 

in a medium with obstacles, correlated at large separations x according to a power law x~ a were analyzed. Integer 
values of the parameter a describe extended defects in the form of lines or planes of random orientation. The obtained 
results qualitatively indicate an increase of shape asymmetry due to the presence of long-range correlated disorder, 
similarly to the case studied in the present paper. In a forthcoming study, we are planning to confirm these results 
numerically and to obtain quantitative estimates of polymer characteristics in the presence of correlated extended 
defects. 



V.B. is grateful for support through the Sachsische DFG-Forschergruppe FOR877 and the National Academy of 
Sciences of Ukraine Committee for young scientists. 
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